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Abstract 

This paper considers a kinetic equation for an unbounded two-dimensional single species point 
vortex system. No axisymmetric flow is assumed. Using the kinetic theory based on the Klimon- 
tovich formalism, we derive a collision term consisting of a diffusion and a drift term, whose 
structure is similar to the Fokker-Planck equation. The collision term exhibits several important 
properties: (a) it includes a nonlocal effect; (b) it conserves the mean field energy; (c) it satisfies 
the H theorem; (d) its effect vanishes in each local equilibrium region with the same temperature. 
When the system reaches a global equilibrium state, the collision term completely converge to zero 
all over the system. 

PACS numbers: 47.32.C-, 47.27.tb, 05.10.Gg, 05.20.Dd, 52.27.Jt, 47.27.-i 



1 



I. INTRODUCTION 



Two-dimensional (2D) microscopic point vortex is a formal solution of the 2D inviscid 
microscopic Euler equation, 

^c2;(r, t) + V ■ (it(r, t)u;{r, t)) = (1) 

where u{r,t) and u{r,t) are the microscopic vorticity and the microscopic velocity, respec- 
tively. Equation ([T]) is formally identical with the macroscopic Euler equation 

^u;(r, t) + V ■ {u{r, t)u{r, t)) = (2) 

where a;(r, t) = {u}{r, t)) and u{r, t) = {u{r, t)) are the macroscopic vorticity and the macro- 
scopic velocity, respectively. Operator (■) is an averaging operator. The point vortex system 
has been successfully applied to the study of 2D turbulence [l, 2|. In the landmark paper 
published in 1949, Onsager proposed an application of statistical mechanics to 2D point 
vortex system, in which he sketched a possible explanation for the formation of large-scale, 
long-lived vortex structures in turbulent flows . Equilibrium distribution at negative 
temperature described by the sinh-Poisson equation is found by Joyce and Montgomery j^. 
Since then, large research effort has been devoted to understand the negative temperature 



state, both theoretically and numerically |6l-ll6j . On the other hand, it has been pointed 



out that a decaying, 2D Navier-Stokes turbulence reaches an equilibrium state described by 



sinh-Poisson equation 17H19| 



Here, a question arises. A distribution of the point vortices oj{r, 0) is given at t = 0. A 
time-evolved distribution uj{r,T) at certain time T is obtained by solving the microscopic 
Euler equation ([1]). On the other hand, a macroscopic vorticity field w(r,0) at t = is 
obtained by a space average of u{r, 0), namely uj{r, 0) = {u}{r, 0)). Of course, a time-evolved 
macroscopic vorticity field u{r,T) is obtained by solving the macroscopic Euler equation (l2l). 
Is the space-averaged point vortex solution {uj{r,T)) the same as the macroscopic vorticity 
filed u{r,T)7 We believe that the answer is "no" and the evolving equation for uj{r,T) 
which is exactly equal to {oj{r, T)) should be written as 

^^^ + V-Hr,tMr,t)) = C (3) 

where C is a collision term. In the following, we restrict our discussion to determining 
explicit formula of C, i.e., a kinetic theory with the Klimontovich formalism j2o| . 



In the context of relaxation processes of the point vortex system towards a statistical 
equihbrium state, diffusive effect due to the discrete distribution has attracted a lot of 
attention. The Fokker-Planck type collision term consisting of the diffusion term and the 
drift (friction) term was discussed by Dubin and O'Neil, in which an axially symmetric 



velocity profile is assumed and a collective effect is formally taken into account |2ll . |22 |. 
This result yields zero diffusion flux for single species point vortices. On the other hand, 
Chavanis obtained the kinetic equation with a collective effect including the implicit term 



23 



24j . By dropping the collective effect, the explicit formula is obtained. However, the 



result has a problem that a monotonic angular velocity profile yields zero coUisional effect 
and the system does not reach a thermal equilibrium state. 



In this paper, we present a kinetic equation for the unbounded, 2D single species point 
vortex system. To derive an explicit formula of the collision term, Klimontovich formalism is 
used. The obtained collision term has following good properties that resolve the above issues. 
The collision term conserves the mean field energy. During a relaxation process towards the 
global equilibrium, system reaches a local equilibrium state first. In the local equilibrium 
state, the relation uj\cq = uj[ip] is satisfied in small regions where w [■?/'] is a functional of the 
stream function ip and in the regions, the second term V ■ {uu) of the Euler equation ([3]) 
vanishes. Then time evolution of the system is dominated by the collision term. However, 
the magnitude of the collision term is small compared to that of V ■ (uu), and the speed 
of the relaxation slows down. When system reaches a global equilibrium state described by 
the Boltzmann distribution Ugeq = ojq exp{— f3Qipgeq) , the collision term completely converge 
to zero all over the system and the Einstein relation is obtained The obtained collision 
term satisfies the H theorem which guarantees that the system relaxes to a global equilibrium 
state. 



The organization of this paper is as follows. In Sec. II, the point vortex system and 
the Klimontovich formalism are briefiy introduced. In Sec. HI, we demonstrate explicit 
formulae for the diffusion and the drift terms as intermediate results. In Sec. IV, a detailed 
calculation of the diffusion term is shown. As the similar calculation can be applied to the 
drift term, details for the drift term are omitted. Finally in Sec. V, three good properties 
of the collision term are demonstrated. 
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II. POINT VORTEX SYSTEM 



Consider a 2D system consisting of positive point vortices. The circulation of each 
point vortex is given by a positive constant Q. 

N 



J2n6{r-n), (4) 



where u = u}{r, t) is the z-component of the microscopic vorticity on the x—y plane, and 6{r) 
is the Dirac delta function in two dimensions. The microscopic variables in the microscopic 
equation are identified by For brevity, we shall omit the t and r dependences, if there is no 
ambiguity. Vector = ri{t) is the position vector of the i-th point vortex. The discretized 
vorticity (jlj) is a formal solution of the microscopic Euler equation ([1]). Other microscopic 
variables are defined by 

ii = u{r,t) = —2 X Vip, (5) 
^ = ^Pir,t) = J2n,Gir-r,), (6) 

i 

G(r) = -i-ln|r|, (7) 

where it and are the velocity field and the stream function in the 2D plane, z is the unit 
vector in the z-direction, and G{r) is the 2D Green function for the Laplacian operator 
with an infinite domain. As a solution of a macroscopic fluid equation should be given by 
a smooth function, the singular solution (jlj) should be regarded as not a solution of the 
macroscopic equation but one of the microscopic equation. Thus we call the equation that 
has the microscopic point vortex solution (jlj), "microscopic" Euler equation. 

Here let us briefly introduce the Klimontovich formalism in plasma physics {20I. It is a 
concise method to derive the Fokker-Planck equation for a macroscopic phase space density 
/ 

^ + vVrf + -{E + « X ■ V./ = V. ■ (D ■ V./ + Af) , (8) 

ot m 

from the Klimontovich equation for a microscopic phase space density /. As the dynamics of 
plasmas is usually dominated by not a collision but a collective behavior due to long-range 
interactions, collision term can be neglected and it yields the simplest form of the kinetic 
equation called the Vlasov equation: 

^ + vVrf + -{E + vxB)-V.f = 0. (9) 
CM m 



We assume that the same hierarchy exists in the 2D fluid equation. The most microscopic 
equation is the microscopic Euler equation ([1]), which have the discrete particle solution (jl]). 
Dividing the microscopic variables into a macroscopic and a fluctuation part, and taking 
ensemble average may yield a macroscopic fluid equation with a collisional effect like the 
Fokker-Planck equation. Ignoring the collision term in the above macroscopic equation, we 
obtain the inviscid fluid equation, namely the macroscopic Euler equation. 

Starting equation is the microscopic Euler equation ([1]). Inserting the following expres- 
sions into Eq. ([1]), 

Co = uj + Su, (10) 
ii = u + 6u, (11) 

and taking the ensemble average, Eq. ([2]) is rewritten as the following macroscopic equation 
with the collision term C = C{r,t) 

^u; + V ■ (uu) = C, (12) 

C = -V-r(r,t), (13) 
r = {6u6uj) 

= - J dr'F{r - r') {6u'6u), (14) 
F{r) = zx VG'(r), (15) 

(16) 

where T denotes a diffusion flux. We note 6u' for 6u{r',t). Similarly, we shall note u' for 
uj{r',t). The following relation is utilized. 

(5w = - y dr'F{r - r')6u' (17) 

In the next section, we will analytically assess the collision term C. 

III. EVALUATION OF COLLISION TERM 

We expect that the collision term C appearing in Eq. f[T^ for the point vortex system has 
two terms, a diffusion term proportional to Vu; and a drift term proportional to u, namely 

T = -D - Vuj + Vuj (18) 



where D = D{r,t) is a diffusion tensor, and V = V{r,t) is a drift velocity. To evaluate D 
and V explicitly, we introduce a small parameter e. We consider a point vortex system with 
large N keeping the total circulation NQ constant. The magnitude of e is of the order of 
either or Q. Orders of the other quantities are: 



C^O(e^), D^O(e), V ^ O(e^), 
uj ^ V'^ifj ^ 0(e°), Vio ^ 0(e), 

VVw ^ 0(e). (19) 



The expansion parameter e is the same as one introduced by Chavanis in Refs. [23|, |2J] and 
the references therein. Expressing D and V in the form of a perturbation expansion and 
gathering the first order terms for D and the second order terms for V, an analytical formula 
for the collision term C will be obtained. 

To rewrite the collision term in Eq. f lT2|) according to the above prospect, we introduce 
a linearized equation obtained by inserting Eqs. ( ITOj) and (ITT!) into Eq. ([1]) and assembling 
the first-order fluctuation terms: 

d 

—Su + V ■ {u6u) = -6u ■ Vu. (20) 

As the macroscopic quantities u appearing in the second term in the left-hand side and 
Vu in the right-hand side are supposed to be constant in the time scale of the microscopic 
fluctuation, Eq. (I20l) can be integrated: 

5u = — drSuir — u{t ~ t),t) ■ Vco 



+5uir-u{t-to),to), (21) 
6u}' = - f dT5u{r' - u'it - r), r) ■ V'o;' 

+5co{r' -u'{t-to),to), (22) 

where V'u' = Vr'u{r'). We call this approximation "straight-line approximation" where 
the trajectory of the point vortex is straight. The value of to is chosen to satisfy t — ^ 
where tc is a correlation time of the fluctuation. Substituting Eqs. fl2T]) and fl22|) into the 
correlation term in Eq. (fT^ . we obtain 

{6u'6u) 



/ dT6u{r' -u\t-T),T) -Vu' 

Jto 

X ^— ^ dr5u{r — u{t — t),t) ■ ViS^ 

^{ir Ik ~ '^'^'^ - ^) ■ '^'^') - ""(^ - ^o), i^o 

+ (5uj{r' - u'{t - to), to) dT6u{r - u{t -t),t) -Vu 



+ {6uj{r' - u\t - to),to)6uj{r - u{t - to), to)) (23) 

— dr {6u{r — u{t — t), t)5u') ■ Vw 

Jto 

- f dr {5u{r' - u'{t - r), t)5uj) ■ W + 0{e^) (24) 

Jto 

^* dr J dr"F{r - u{t - r) - r") ■ Vu {6uj{r", t)6uj') 

+ I'drJ dr"F{r' - u'{t - r) - r") ■ Vu' {6uj{r", t)6uj) . (25) 



When obtaining formula ( 124|) . we assume that the first term in formula ( 123|) is negligible as 
it has two nablas. We drop the last term as it should have a factor of l/(t — to) and we focus 
on t — to ^ case. The time is shifted from to to t using the straight-line approximation. 
When rewriting formula flMj) as fl25|) . Eq. f|T7|) is used. 

Inserting Eq. ( l25l) into Eq. ( fT4l) . the following intermediate results are obtained: 



D ■ Vw 

ft 



f dr f dr' f dr"F{r - r')F{r - u{t - r) - r") ■ Vw 

Jto J J 

X {5u{r",T)6uj') , (26) 
Vio 

- fdrj dr' J dr"F{r - r')F{r' - u'{t - r) - r") ■ VJ 

X {5u{r",T)5uj) . (27) 



It should be noted that the diffusion term can be expressed in the following modified Kubo 
formula: 

D= dT{6u{r,t)6u{r - u{t - t),t)). (28) 

Jto 



IV. EVALUATION OF DIFFUSION AND DRIFT TERMS 



As the expression of the diffusion term (l26l) is very similar to that of the drift term (127|) . 
the detailed derivation for diffusion term only is shown: 

{6uj{r'\ t)6uj') 
= {[u{r'\r)-uj{r'\r)] [uj' - u']) 
= {Cj{r'\T)Cj') - uj{r\T)J 

INN \ 
\i=l jjti I 

-u:{r\T)u'. (29) 

The first term in the last result in Eq. fl29|) corresponds the case of i = j, and the second 
term corresponds the case oi i ^ j. 

For the i = j case, the formula is rewritten as 

= (j2 ^'Kr" - nir) -r' + n{t))6{r' - r,(t))^ 

N 

= J2^' {S{r" - r.{r) - r' + r,{t))5{r' - r,(t))) . (30) 
1=1 

Here we introduce a stochastic process to evaluate rj(t) — rj(r): 

r,(t) - n(r) = J' u(r(r'), T')dT' + ^ 

^ u'it -T) + i. (31) 

The first term in Eq. flSTl) represents the straight-line approximation and the second term 
represents a Brownian motion. The stochastic process represented by (■)^ includes all the 
possible motion to reach position at time t. Then, Eq. ( l30i) can be rewritten as 

N 

Y: {Sir" - r,(r) - r' + r,(t))5(r' - r,(t))) 

i=l 
N 

= J2^' {5{r" - r' + u\t - r) + 0)^ (5(r' - r,(t))) 

i=l 

= {5{r" -r' + u\t-T)+i))^^uj'. (32) 
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For the % ^ j case, we introduce an approximation that correlation between the particles 
can be neglected 

N N 

EE^' {S{r" - r,{r))S{r' - r,m 

i j^i 
N N 

^EE^' {Kr"-n{r))) (5(r'-r,(t))). 
Also we assume the followings: 

N 



i 

= NQ {6ir" - r,(r))) . 



(33) 



(34) 



Inserting Eq. into Eq. (1331) . we obtain 



N N 



i j^i 



/ N \ ( ^ ^ 

[Y.^{5{r"-r,{r)))\ x [Y.^ {5{r' - r,{t))) 
\i=i I . 



NQ 



NQ 



= co{r",TW--uj{r",T)u'. 
Combining the results of i = j and i j cases, we rewrite Eq. (129|) as 



(35) 



{6uj{r",T)5iu{r',t)) 

Q{6{r" -r' + u'{t-T) +^))^u' 



(36) 



The two terms in the right hand side of Eq. (!36|) are the same order terms as we request 
the total circulation NQ constant. To proceed with the evaluation, a conservation law is 
introduced 

J dr' {6uj{r",T)6u') =0. (37) 
Inserting Eq. fl36|) into Eq. fl371) . we obtain 

dr' [Q {6{r" -r' + u'{t - r) + ^))^ u' 
^ (II \ I 
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This equation yields 



Q J dr' {5{r" - r' + u'{t - r) + $,))^ u' 

-^uj{r",r)jdr'uj' 

Q J dr' {6{r" - r' + u\t - r) + i))^ J 

0. (38) 



I dc^ {5{r" -q' + u{q') {t - r) + i))^ u{q' , t) (39) 



where dr' is replaced by dq' to avoid ambiguity. This equation enables that all the quantities 
at r is converted by ones at t. Inserting Eqs. (|36|1 and (l39l) into Eq. (I26l) . we obtain 



D ■ Vw 

rt 



Q f dr f dr' f dr"F{r - r')F{r - u{t - r) - r") ■ Vu 

Jto J J 

x{S{r"-r' + u'{t-T) + ^))^iu' 

~Jl It I I ^^"■^^'^ ~ ^')^(^ - u{t - r) - r") ■ Vu 

xu' J dq' {5ir" -q' + u{q'){t - r) + ^))^ io{q' , t). (40) 



We proceed with the evaluation of the second term in Eq. fHOl) : 

-— r dr f dr' f dr"F{r - r')F{r - u{t - r) - r") ■ Vu 
xu' J dq' {5{r" -q' + u{q'){t - r) + ^))^ uj{q' , t) 

= -l| dr'F{r-r')j'^dr J dq' {F{r - q' - {u - u{q')){t - r) + ^))^ ■ Voo 

xuj'uj{q',t) (41) 

1 /■ , /,-,/ , f dk z X ik ^ 

= -— / dr'F(r - r') / dr / - — - ,^ ■ Vu 

NJ ^ ^ Jto J (27r)2 |fc|2 

X J dq' ex.p[ik ■ {r — q' — [u — u{q')){t — r))] (exp(2fc ■ $))^ 

xuj'uj(q',t). (42) 

To rewrite formula (HTl) as ( H2l) . we use the Fourier transformation: 

i^(r--g'-(n-ii(g'))(i-r)) 
= (2^ / dk^^expitk-ir-q'-iu-uiq'))it-r))). (43) 
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The term (exp(ifc ■ ^))^ represents a Brownian motion of the point vortices with diffusion 
tensor D and is evaluated by the cumulant expansion: 

k-D-k 



{exp{ik ■ = exp 



-it-r) 



= exp{-u{t - r)) (44) 
where u is a small positive parameter. Inserting the following formula into Eq. ( H2l) 

/ dr exp[—ik ■ {u — u{q') — v){t — r)] 



n5{k ■ {u - u{q'))) 



ik ■ {u — u{q')) 
\k-{u-u{q'))\^ + iy^' 



(45) 



we obtain 



~ I dr'F(r-r') t dr 



dk z X ik 



■ Vu 



X J dq' exp [ik ■ {r — q' — [u — u{q')){t — r))] (exp(ifc ■ ^))^ 
xuj'uj{q', t) 



1 1 dr'F{r - r') j 



dk z X ik 



N 
X I dq' 



TT 



5{k- {u-u{q')) 



ik ■ [u — u{q')) 
\k ■ {u - u{q'))\^ + 



exp{ik ■ (r — q')) 



xu'u{q', t). 



(46) 



We substitute r + q" for q' and expand u{q') and uj{q') in the form of Taylor series and 
retain the zero-th order terms only: 



u{q') = u{r) + q" ■ Vu{r) + 0(e), 
tuiq') = u{r + q") = uir) + Oie). 

Inserting Eqs. (14 7p and (1481) into Eq. (H^ . we finally obtain 



(47) 
(48) 



-— / dr'F(r - r') 



dk z X ik 



(27r)2 |fc|2 



X / dq' 



7r6{k ■ {u — u{q')) 



xu'u{q', t) 
1 



X / dq" 



dr'Fir - r') 



ik ■ (n — u{q')) 
\k ■ {u - u{q'))\^ + 

dk z X ik „ 



exp{ik ■ (r — q')) 



n5{-k ■ {q" ■ V)u) 



ik ■ {q" ■ V)u 
(q"- V)n|2 + i/2 



exp(— ifc ■ q") 



xu u. 



(49) 
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It is found that Eq. (jUj) substituting k = —k and q" = —q" changes its sign. Thus 
it is concluded that the integral equals zero, i.e. the second term in Eq. fHOj) has zero 
contribution and only the first term remains. The obtained formula for the diffusion term 
is as follows: 



-D ■ Vu 



n J dr'F{r - r') j exp{ik ■ (r - r'))^^4^ ■ u'Vu 



k 



7i5{k ■ {u — u')) 



ik ■ iu — u') 



(50) 



\k ■ {u - u')\^ + e\ 

The similar calculation can be adapted for the drift term. For this case, the following 
conservation law is used: 

J dr {6io{r", r)5w(r, t)) = 0. (51) 
The whole result including both diffusion and drift terms are given by 



-D ■ Vw + Vu 
Q I dr' ' '^^ 



dk' 



(27r)2y (27r) 



exp(i(fc + k') ■ [r — r')) 



nS{k ■ (u — u')) 



ik ■ {u — u') 
\k ■ {u — u')p + 



z X xk z X xk , , , ,, 

X — ^-t:^ r^-:^ ■ {u'Vuj - uV u') 



\k 



i\2 



\k\ 



(52) 



where we have used Eq. 



V. SPACE-AVERAGED COLLISION TERM 

Equation (152|) includes the oscillatory term exp(i(fc+fc')-(r— r')). To reveal characteristics 
of the obtained collision term, we need to calculate the space average of the collision term to 
drop the high-frequency component. Space average is calculated over the small rectangular 
area A with sides both 2L located at r. The space average of the diffusion flux T given by 
Eq. ( 15^ is defined by 

We assume that the macroscopic variables such as u and u may be constant inside A(r) 
and only the term exp(z(A; + k') ■ [r — r')) should be space-averaged: 

(exp(^(fc + k') ■ {r - r')))s 

12 



(2L) 



^ dx" dy"exp(i(k + k') ■ r") exp(-i(k + A;') ■ r') 

LV J-L J-L 



-) 5(fc + fc')exp(-z(fc + A;') - r) 



where r" = Thus, space-averaged diffusion ffux is given by 



(54) 



(2vrj 

z k z k . , , 



In Eq. fl55|) . we omit the imaginary part as the collision term consists of only the real part. 
Further integration over k in Eq. (155|) can be performed. The integral concerning k is as 
follows: 

, [z X k){z X k) 



J dk6{k ■ {u - u'))- 



\k' 



(56) 



Dividing k into the parallel and the perpendicular components and inserting them into Eq. 

m, 

k = k\\ri^\ + k±h±, 



u — u 



nil 



\u — u'\ 
n_|_ = I; X nil, 



(57) 



we obtain 



(z X k){z X k) 



J dk6{k . {u - u')) 
dkii J dk_L5{k\\\u — u'\) 



X 



[z X {k^\ri\\ + kj_hj_)][z x {k\\n\^ + kj_nj_)] 



2 12 



dkj 
dk 



1 



\kl + k 
[z X k±h±)(z X k±hj_) 



\u — u'\ k'\_ 

1 1 u — u' 1 u — u' 



\u = w 



i'\kl 



U — U' \ K, \U — W 



(u — u')(u — u') 



\u — u 



'13 



{u-u'){u-u') 2 



\u — tt'P 



kn 



(58) 
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where parameter fcmin is introduced to regularize a singularity and determined by the largest 
wave length that does not exceed a system size, namely fcmin = 27r/i? where i? is a charac- 
teristic system size determined by an initial distribution of the vortices. 
Finally, we obtain the following formulae for the diffusion and drift: 

T,ir) = -Ds{r)-Vco + Vs{r)u, (59) 
n /7r\2 1 f ,(u - u')(u - u')uj' , 
■^'-(2^(1) k-.h'- lu-^f ■ 

(27r)3 \lJ k^i^J \u-u'\^ ^ ' 



In Eqs. fl60|) and fl6T|) . two unknown parameters L and /cmin remain. We assume that L = gR 
where 5^ is a size factor (^f ^ 1). If we set (7 = l/47r, Eq. (159!) is rewritten as 

^ , , Q r , ,(u — u')(u — u') 



R J \u — u'\^ 

u'Wcu -tuVu'). (62) 



A. Collision term in local and global equilibrium states 



At first, let us examine if the coUisional effect (16211 disappears in a local equilibrium state. 
We rewrite Eq. ( l62l) into a symbolic form: 

T,ir) = -^ldr'j[u,t,oj',^'] (63) 

where 7 is a functional of u, ip, u' and ip' . Consider a state where temperature is locally 
uniform in each small region in the system. We call this state the local equilibrium state in 
which the local equilibrium condition is satisfied: 

Wleq = Wo exp(-/3il?/'ieq). (64) 

Inserting Eq. ( I64l) into 7 in Eq. ( I62l) . we find that 



7[u;ieq,^leq;U;i'eq,^leq] 

(''^leq ~ '^leq) / , . , , ^ , ^ 

, ig l^leq - ^^leq) " (^IcqV^leq - WjcqV Wj^J 

I "Icq "Icq I 

-/JfiWleqWi'eqi^^^i^^— ^^(Wlcq " w(cq) " (V^lcq - VVl'eq) 
[""leq ""leql 



(65) 
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where uicq = —z x V^^icq is used. As uieq — w(eq is perpendicular to VV'icq — VVi'eq; 7 is equal 
to zero and this result indicates a detailed balance is achieved. When the system reaches a 
global thermal equilibrium state jsl 

CJeq = exp(-/3n?/'eq) , (66) 

we obtain 

Vc^eq = ^eq— 

eq 

= -/3fiu;^q(VV:q - V^eq + V^^eq)- (67) 

As (Ueq — u'^q) ' (V^eq ~ V'i/'eq) = 0, the drift tcrm in Eq. (!62l) is rewritten as 

V,,eq = -/3^^D,,eq " V^eq (68) 



which is a counterpart of the Einstein relation 23|. 



B. Energy-conservative property of collision term 

Time derivative of the total mean field energy E is given by 



^='-fdrfdr'Gir-r') (^4u + u'^] 
dt 2J J ^ ' \dt dt j 



rfr^— , (69) 



where 



= ^Jdr J dr'G{r - r')uj'u. (70) 

Note that the mean field energy E is different from the system energy ^ of the point vortex 
system 

M^ = -^T.ll^^^Mr^-r^. (71) 
Inserting the space-averaged equation of motion 

^ + V ■ {uu) = -V ■ r, (72) 
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into Eq. (l69i) . we obtain 
dE 



J drip (- V ■ (uu) - V ■ r,) 
J drVip ■UU + J drVi) ■ 
J drVip ■ r. 



^ /■ , /■ , ,„ (u — u')(u — u') , . „, , ^ 

KJ J lu — u'l-^ 

By permuting the dummy variables r and r' in Eq. fl73l) and taking the half-sum of the 

resulting expressions, we obtain 

X (it - u') ■ {uVu - uVu') 



\u — u'P 



= 0. (74) 
We conclude that the obtained collision term conserves the total mean field energy. 



C. H theorem 



The entropy function S is defined by using the H function: 

S = -ksH, (75) 

f ^ u ^ u 
±/ = / dr— In — + const. 

= — J dru Into ~ N Inil + const.. (76) 



The time derivative of the H function is given by 

= dri — (lnw + 1 
dt nJ \ dt^ ' 



rfr(-V - (uw) - V-r,)(lnw + l) 
^ / druoj ■ V In + / drFg ■ V In w 

I- f dr{V ■u)uj + ^ f drT.-Vlnu 

\l J il J 



n 

^1 drT.-Vlnu. (77) 



Inserting Eq. ( I62l) into Eq. ( ITTl) . we obtain 



^ Urfdr'^'' (n-n')(--i.' 



dt RJ J uj \u — u'\^ 

■{u'Vuj -uVuj') (78) 
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By permuting the dummy variables r and r' in Eq. (178|) and taking the half-sum of the 
resulting expressions, we obtain 



dH 1 /",/■,, 1 , ^/ /N u-u 

— = / dr / dr' -(coVtu - wW) ■ -. — 



dt 2RJ J uu''" "' \u-u'\^ 
X (it - u') ■ (w'Vw - ujVu') 

iu-u') ■ {u'Vu -uVu')? 



-— I dr I dr'— 
2RJ J uouj' 



\u — u 



/|3 



< 



(79) 



The integrand of Eq. fl79|) is positive or equal to zero, and dH / dt is negative or equal to zero. 
It is concluded that the entropy function S ( |75l) is the monotonically increasing function. 



VI. DISCUSSION 

We have demonstrated the simple and explicit formula of the Fokker-Planck type collision 
term for nonaxisymmetric point vortex profile without the collective effect. It should be 
noted that it does not include the implicit term like Eq. (26) in Ref. 2J]. The collision term 
exhibits several important properties: (a) it includes the nonlocal, long-range interaction; 
(b) it conserves the mean field energy; (c) it satisfiies the H theorem; (d) its effect vanishes 
in each local equilibrium region with the same temperature. When the system reaches a 
global equilibrium state, the collision term completely converge to zero all over the system. 

There are several outstanding issues remaining. First, the final formulae (160!) and (!6T|) 
includes unknown parameters kmin and L. Second, the integrals in Eqs. fl60l) for D and (|6T|) 
for V contain the divergent integrand, although the combined term (r)^ = — D ■ Vw + Vuj 
is regularized. A more rigorous justification will be needed for fixing the above two issues. 
Finnaly, although our result can be adopted for an axisymmetric case, it does not yield the 
same results presented by Dubin, O'Neil 
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22| and Chavanis [23|, |2J] . A reason remains 



unknown yet. 
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